Effective Hamiltonian for Gai_ x Mn x As in the Dilute Limit 
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We derive an effective Hamiltonian for Gai-^Mn^As in the dilute limit, where Gai-^Mn^As can 
be described in terms of spin F — 3/2 polarons hopping between the Mn sites and coupled to the 
local Mn spins. We determine the parameters of our model from microscopic calculations. Our 
approach treats the large Coulomb interaction in a non-perturbative way, captures the effects of 
spin-orbit coupling and disorder, and is appropriate for other p-doped magnetic semiconductors. 
Our model applies to uncompensated Mn concentrations up to x ~ 0.03. 

PACS numbers: 75.30.Ds, 75.40.Gb, 75.50.Dd 
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Since their discovery^, dilute III-V magnetic semicon- 
ductors with high Curie temperatures have become the 
subject of very intense research Because the magnetic 
ions (usually Mn) responsible for the ferromagnetism are 
dissolved into the semiconductor itself, these materials 
could provide a unique opportunity to integrate ferro- 
magnetic elements into larger, non-magnetic, semicon- 
ducting devices. 

In this Letter we focus on one of the most studied mag- 
netic semiconductors, Gai-^Mn^As, though most of our 
calculations carry over to other p-doped III-V magnetic 
semiconductors |3j. In Gai-^Mn^As substitutional Mn 2+ 
play a fundamental role: They provide local spin S — 5/2 
moments, and they dope holes into the lattice^]. Since 
the Mn 2+ ions are negatively charged compared to Ga 3+ , 
in the very dilute limit they bind these holes, forming an 
acceptor level with a binding energy Eb ~ 112meV 0. 
As the Mn concentration increases, these acceptor states 
start to overlap and form an impurity band, which for 
even larger Mn concentrations merges with the valence 
band. Though the actual concentration at which the im- 
purity band disappears is not known, according to op- 
tical conductivity measurements^, this impurity band 
seems to persist at least up to nominal Mn concentra- 
tions of about x w 0.05. ARPES dataQ and the fact 
that even "metallic" samples feature a resistivity upturn 
at low temperatureQ suggest that for smaller concentra- 
tions (and maybe even for relatively large nominal con- 
centrations) one may be able to describe Gai-^Mn^As 
in terms of an impurity band ||. 

In Gai-^Mn^As the Coulomb potential created by the 
Mn ions is by far the largest energy scale in the prob- 
lem ||, but spin-orbit coupling in the hole band is also 
quite large compared to the exchange coupling between 
the holes and the Mn spins 0. Fortunately, the large 
Coulomb potential of the Mn ion can be handled non- 
perturbatively. We construct a many-body Hamiltonian 
in this limit that captures spin-orbit effects, treats the 
large Coulomb interaction non-perturbatively, and in- 
corporates the exchange coupling between the local mo- 



ments and the holes. 

The physics of the isolated Mn + + hole system is 
well-understood jj]: In the absence of the Mn 2+ core 
spin, the ground state of the bound hole at the acceptor 
level is four-fold degenerate and well described in terms 
of a F = 3/2 spin. For most purposes, we can restrict 
ourselves to this fourfold degenerate F = 3/2 acceptor 
level in the dilute limit. As also evidenced by infrared 
spectroscopy |^|, the effect of the S — 5/2 Mn core spin 
is well-described by a simple exchange Hamiltonian |2(: 
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with G w 5 meV |4|. 

The presence of other Mn sites has three important 
effects on the F — 3/2 acceptor state at any particular 
Mn site: (i) The Coulomb potential of the neighboring 
Mn 2+ ions will induce a random shift E of the fourfold 
degenerate states, (ii) Because of the large spin-orbit 
coupling in GaAs, the neighboring atoms will also gen- 
erate an anisotropy K and split the fourfold degeneracy 
of the F = 3/2 state into two doubly degenerate states. 
(Hi) Finally, the presence of the neighboring ions will al- 
low these F = 3/2 spin objects to hop between the Mn 
sites. However, this hopping t will not conserve the spin 
F because of the large spin-orbit coupling. Thus, in the 
dilute limit Gai-^Mn^As should be described by the fol- 
lowing simple Hamiltonian: 

H = E^c^+E^^r+^n^ 
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where c\ v creates a hole at the acceptor level \F = 
3/2, F z = v) at position i. As is clear from the argu- 
ments above, Eq. (0) is very general and appropriate for 
describing other p-doped III-V and II- IV semiconductors 
as well [3j . Hole-hole interactions can be incorporated in 
Eq. J5J). In the metallic phase it is presumably a good 
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approximation to include only an on-site repulsion (dis- 
cussed later) that eliminates double occupancy of the ac- 
ceptor levels. In the localized phase, however, one may 
have to consider long-ranged hole-hole interactions. 

To estimate the various parameters in Eq. we stud- 
ied the structure of the single impurity (Mn) and two- 
impurity (MJI2) bound hole states using the so-called 
spherical approximation 10J. The top of the valence band 
in Ga i -y Mn^As can be described in terms of spin j = 3/2 
holes , whose spin couples strongly to their momenta. 
In the spherical approximation the motion of the holes in 
the Coulomb potential of an Mn ion is described bv[Io|. 
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where 7 « 7.65 is a mass renormalization parameter, m 
is the free electron mass, n « 0.77 is the strength of the 
spherical spin-orbit coupling in the j = 3/2 band |l0| . 
e w 10 is the dielectric constant of GaAs, and V cc is the 
so-called central cell correction ^| . The spin-orbit term 
in Eq. J3J couples the momentum tensor of the holes 
Pa/3 = PaP/3 — S a p p 2 /3 to their quadrupolar momentum, 

7 Q/ 3 = (j a jp +J>ia)/ 2 - <W j(j + l)/3. 

The bound states of H (without the central cell cor- 
rection) have been studied in the seminal paper |Toj. 
Due to the spherical symmetry, the total momentum, 
F = L + j, is a conserved quantity, where L is the orbital 
angular momentum. The ground state of Hq is a four- 
fold degenerate F = 3/2 multiplet that contains a sub- 
stantial d-wave contribution for GaAs due to the strong 
spin-orbit coupling. In Fig. ^ we illustrate the impor- 
tance of this d-wave component by presenting the spatial 
dependence of the hole polarization, j(r), for the state 
\F = 3/2, F z = 3/2) which we calculated directly from 
the Baldereschi-Lipari wave functions with the central 
cell correction. 

Having computed the single Mn hole states, we carried 
out a variational calculation to construct the molecular 
orbitals for a pair of Mn ions 0, Q] . Since the exchange 
interaction with the Mn core spins is much smaller than 
the binding energy of the holes, we neglected the effect 
of G on the their wave functions in these calculations. 
For the case where both the Mn-Mn bond and the quan- 
tization axis of F are parallel to the z-axis, F z is con- 
served and the spectrum of the lowest lying states of the 
molecule can be fully characterized by: 
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By time reversal symmetry, the hopping parameters sat- 
isfy t 3 / 2 = £-3/2 an d = t-i/2- All parameters de- 
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FIG. 1: Polarization of a bound hole in the state \F = 
3/2, F z = 3/2} in Gai_ x Mn x As around a Mn ion (dark arrow 
pointing downwards represents the Mn S=5/2 spin). Only 
the direction of the polarization is indicated. The magnitude 
falls off on a scale ~ 10 A. 
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FIG. 2: Parameters of the two-impurity Hamiltonian Eq. Q 
obtained from the variational study of two Mn ions. The 
arrow indicate the typical Mn-Mn distance, dt yP , for x = 0.01 
Mn concentration. 



pend only on the distance R between the two Mn sites 
(see Fig. |2J). The most obvious effect of the spin-orbit 
coupling is that the hoppings £3/2 and t\/2 substantially 
differ from each-other; holes that have their spin aligned 
with the Mn-Mn bond are more mobile. As indicated by 
the arrow, at the typical Mn-Mn distance for x = 0.01, 
K and t\ii can be entirely neglected compared to E and 
t 3 /2- Therefore, in many cases it is enough to keep only 
the latter two terms in the effective Hamiltonian. 

Having determined the effective Hamiltonian for a 
pair of Mn ions, we can use it to estimate the pa- 
rameters in Eq. (J2J. Rotating the z-axis along the 
bond direction Hij connecting sites i and j, we obtain 
U, = P(n y )t(i? y )I?t(^.) ; Kl = IZ^K(R tJ ){(n^ ■ 



3 



Ff - |), and Ei = §£ j¥i -E(i2«)- Here D(fty) is 
a spin 3/2 rotation matrix, and t(R) denotes the di- 
agonal matrix diag (t 3 / 2 (R), t 1/2 (R), t 1/2 (R), t 3/2 (R)), 
and Rij denotes the distance between sites % and j. 

So far we have neglected the interaction between holes. 
In the localized phase, however, this interaction may play 
an important role. In general, the Coulomb interaction 
between holes on different Mn sites has a very compli- 
cated form |3|, though for large separations it simplifies 
considerably. Fortunately, the dominant interaction is 
the on site hole-hole interaction. Within the spherical 
approximation this interaction can be expressed as: 

i i 

where Ni = Z)„ ct F t = q^i^Ci,^ and : ■■• : 
denotes normal ordering. We estimated Un and Up in 
Eq. Jfjjl by evaluating exchange integrals: Un — 2600K 
and U F = -51 K. 

Eqs. J21 and (J5J, together with the microscopic param- 
eters of Fig. |3 constitute our central results. They pro- 
vide a well controlled theoretical framework that captures 
the most important aspects of dilute magnetic semicon- 
ductors such as the localization phase transition, ran- 
dom anisotropy, disorder effects, and frustrated ferro- 
magnetism. Postponing much of our detailed analysis 
to a longer publication |3J, here we only demonstrate the 
power of this model on a few examples. 

To obtain a better understanding of the model we first 
computed the ground state of four Mn atoms at a sepa- 
ration of 15 A due to the interaction mediated by a sin- 
gle hole on the cluster. We treated the Mn spins clas- 
sically and used the simple mean field approximation of 
Ref. |8| • We considered only configurations where the Mn 
ions were positioned on a slightly distorted tetrahedron 
with three edges of length a = 15 A and three edges of 
length b (see Fig. EJ). In all cases, in the ground state, 
the Mn spins are relatively collinear apart from a slight 
tilt of 5 — 10°. However, the spatial position of the Mn 
ions generates a strong anisotropy. Thus, the energy de- 
pends strongly on the directional orientation of the net 
spin relative to the underlying lattice. To demonstrate 
this we calculated the ground state energy as a function 
of the Mn spin direction (assuming full alignment) . For 
a perfectly regular tetrahedron this anisotropy is rather 
small, less than 0.5K/Mn. However, the anisotropy in- 
creases with the ratio b/a, and for b/a = 1.2 it can be 
as large as 20K/Mn. In other words, random positions 
of the Mn ions induce a random anisotropy term that, 
depending on the disorder, is much larger than the bulk 
anisotropy, which is of the order of lK/Mn. Thus disor- 
der and spin-orbit coupling together can induce a large 
random anisotropy energy comparable to Tq- These find- 
ings are in good agreement with earlier results obtained 
in the metallic limit 1151. 
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FIG. 3: Anisotropy induced by the distortion of a regular 
Mn tetrahedron in the presence of a single hole. The Mn- 
Mn distances are a = 15A and b = 15A, 6 = 16. 5A, and b = 
18A, respectively. The distortion generated anisotropy can be 
almost two orders of magnitude larger than the undistorted 
anisotropy, which is of order lK/Mn. 



Finally, we discuss some of the results obtained for a 
Ga!_ x Mn x As of linear sizes L = 10ai at and L = 13ai at 
(with ai a t = 5.65A, the size of the conventional unit cell) 
and active Mn concentration ^active = 0.01, using the 
above-described mean field techniques at zero tempera- 
ture. In the calculations presented below, we have not 
included the effects of Eq. JHJ. (This is partially justified 
post-facto by Fig. 01 which shows the states at the Fermi 
energy are delocalized. If these states were localized one 
would expect the contributions of Eq. JSJ) to be more im- 
portant.) We emphasize that x active can be substantially 
less then the nominal Mn concentration, x, which also 
includes inactive Mn sites and therefore these cal- 
culations may be relevant even for systems with larger 
nominal Mn concentration. The concentration of holes is 
also reduced compared to x due to strong compensation 
effects; we assumed that the number of holes is reduced 
by a factor of / = 0.3 relative to the number of Mn. 

To take into account correlations induced between 
Mn ions during the experimental growth process, we in- 
troduced a screened Coulomb repulsion between the Mn 
ions and let them relax using T=0 Monte Carlo simu- 
lations. For long times the Mn ions form a regular bcc 
lattice with some point defects. The data we present here 
are for intermediate times, where there is still appreciable 
disorder in the system. 

Once we fixed the Mn positions in a given instance, we 
solved the mean field equations derived from self con- 
sistently 0- We used periodic boundary conditions and 
implemented a short distance cutoff in the hopping pa- 
rameters of Eq. J5J which corresponds to about 8 neigh- 
bors for each Mn. The use of this cut-off is justified by 
the observation that our molecular orbital calculations 
are only appropriate for "nearest neighbor" ion pairs, 
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FIG. 4: Top: Computed average hole density of states for 10 
L = 10a samples with x = 0.01 and / = 0.3. We also show 
the density of states of the valence band (dashed line). The 
Fermi energy is « —6500 K. Bottom: The participation ratio 
for L = 10ai a t and L — 13ai at . States in the impurity band 
tails are localized while states in the middle are delocalized; 
states in the side tail above zero energy likely mix with the 
valence band states and are delocalized in reality. 

and in reality, holes cannot hop directly over the first 
"shell" of ions. We started from a configuration with 
fully alligned classical Mn spins, f2j = S/S, and then let 
the system relax to the nearest metastable state. Similar 
to the metallic case 01 > we find a ferromagnetic state 
with a largely reduced magnetization, |(f2j)| w 0.4 for 
L = 10ai a t. We find that this reduction is largely due 
to spin-orbit coupling, and that the cosine of the angle 
9 between the spins and the ground state magnetization 
has a broad distribution similar to the metallic case . 

The density of states is shown in Fig.0] The halfwidth 
of the impurity band is about 0.2 eV at this density, 
which slightly overlaps with the valence band density of 
states. However, comparison with the valence hole den- 
sity of states suggests that at this concentration a well- 
formed impurity band may still be present, and it might 
persist to higher concentrations. Indeed, this scenario 
seems to be supported by many experiments|5j,|6|]. 

The impurity band has a tail of localized states that 
reaches inside the band gap. These states can be identi- 
fied for various system sizes. (See Fig. 0] wherein the par- 
ticipation ratio, PR = K,-CCc« l^ia| 2 ) 2 ]~\ grows with 
system size for delocalized states while the PR remains 
0(1) in the thermodynamic limit for localized states). 

This tail gradually disappears when we introduce cor- 
relations between the Mn ions which tend to form reg- 
ular structures 0]. In agreement with ARPES dataQ, 
we find that the chemical potential lies deep (~ 0.5 eV) 
inside the gap. From the PR data, it appears that the 
chemical potential is in the vicinity of the mobility edge, 
a regime where our model is probably more reliable. 



This raises the interesting possibility that the localiza- 
tion phase transition in Gai^^Mn^As could happen in- 
side the impurity band and that the ferromagnetic phase 
for smaller Mn concentrations is governed by localized 
hole states 0, . 

Though our calculations are based on microscopic 
model calculations, they are only approximate, and more 
realistic ab initio calculations would be needed to give a 
quantitative answer concerning the role of the impurity 
band. Also, though the spherical approximation we used 
is able to reproduce rather well the spectrum of a single 
acceptor, it might overestimate the effect of spin-orbit 
coupling, and also the width of the impurity band. 

In summary, based on microscopic calculations we con- 
structed a many-body Hamiltonian that is appropriate 
for describing Ga^^Mn^As in the dilute limit. We find 
that the hopping of the carriers is strongly correlated 
with their spin. This spin-dependent hopping is cru- 
cial for capturing spin-orbit coupling induced random 
anisotropy terms, the lifetime of the magnon excitations, 
and for capturing the universality class of the localization 
phase transition. Our calculations suggest the presence 
of an impurity band for rr a ctive = 0.01 Mn concentration. 
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